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Abstract 

We show that an unambiguous way of determining the universal hmiting fragmentation region 
is to consider the derivative {d?n/drf') of the pseudo-rapidity distribution per participant pair. In 
addition, we find that the transition region between the fragmentation and the central plateau 
regions exhibits a second kind of universal behavior that is only apparent in d^n/drj^. The y/s de- 
pendence of the height of the central plateau {dn/dri)ri=o and the total charged particle multiplicity 
ntotai critically depend on the behavior of this universal transition curve. Analyzing available RHIC 
data, we show that {dn/dr])j^=Q can be bounded by In^ s and ntotai can be bounded by In^ s. We 
also show that the deuteron-gold data from RHIC has the exactly same features as the gold-gold 
data indicating that these universal behaviors are a feature of the initial state parton-nucleus in- 
teractions and not a consequence of final state interactions. Predictions for LHC energy are also 
given. 
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I. INTRODUCTION 



Recently, PHOBOS and BRAHMS collaborations at RHIC published a set of intriguing 
data. Among them are the striking feature of the limiting fragmentation It is reported 

by PHOBOS that the shifted (by the beam rapidity ^max) pseudo-rapidity distribution per 
participant pair dn/dr] = ((iA^ch/'^^)/(A^part/2) is independent of colliding energy up to 85- 
90 % of the plateau height Hi. 

If taken at face value, this would imply that the height of the plateau and its ^/s de- 
pendence would be almost fully determined by the limiting fragmentation curve and the 
location of the beginning of the plateau. Also, since the total multiplicity is simply the area 
under dn/drj, its y/s dependence would be largely determined by the limiting fragmentation 
curve as well. 

It is not easy to determine the validity of above statements when only the rapidity 
distributions {dn/drj) are compared. In this paper we argue that comparing the slopes 
[d^n/drf) is a much better way to determine various regions of the rapidity distribution. 

A surprising feature is the existence of another 'universal' behavior, which is only appar- 
ent in the slopes. It turns out that d'^n/drf in the transition region between the limiting 
fragmentation and the central plateau also follows a universal curve. This is not an extension 
of the limiting fragmentation curve. To our knowledge, this is the first time the existence of 
this second universal curve is demonstrated. These two universal curves basically determine 
the dn/drj. The energy dependence shows up through the position of the beginning of the 
central plateau. 

The hypothesis of limiting fragmentation has a long history. For hadron-hadron collisions, 
this hypothesis was first put forward by Benecket et.al.|iand also by FeynmanW and 
Hagedorn ^ . This idea was further developed in Refs. ^ . 

Feynman hypothesized that as ^/s oo, the multiplicity spectrum 

lim £;,§^= lim ;^ = /(.„p,) (1) 

y/s^oo d-^p v^-»oo dy d'^pT 

becomes independent of -/s. Here y = (1/2) ln((i?p + pl) /{Ep — pl)) is the rapidity and 
xl = 2^^/1/5 is the longitudinal momentum fraction. If the mass of the particles is light 
compared to the average pt, this expression also equals drihh/ dr] d'^px where rj = (1/2) ln((p+ 
Pl) I {p—pi)) is the pseudo-rapidity. The universal function f^xi^Pr) then totally determines 
the height of the dnhh/drj and the total multiplicity at high energies. 

Note that since f^XL^Px) itself is independent of y^, the height of the plateau 
[duhhl dy)^^Q must also be independent of i/i. This also implies that the total multiplicity 
must behave like ?/max ~ In s where 

1/max = COsh~^{^/s/2mN) ~ \n{y/~s/mN) (2) 
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is the beam rapidity and mjv is the nucleon mass. However, up to = 1800 GeV the 
experimental data does not show that {dn/d'r])rj=o is saturated. Also proton-proton and 
proton-anti-proton data show that the height of the plateau grows like In^ s (See compilation 
by PHOBOS in Refs.QllZl .). The the total multiplicity then must grow like In^ s. This is 
not what one would expect from Eq.ijT]). 

The source of this discrepancy is the fact the strict Feynman-Yang scaling is not perfect 
nor is it supposed to be. The central region (or small xl) is modified by radiation of soft 
partons and and the multiple rescatterings of produced particles. QCD radiative corrections 
should also give rise to the additional scale dependence in f(xr. pt-)\i^. 

However, within the dynamic range where Feynman-Yang scaling approximately holds, 
what should still work is the universality of dn/drj near rj = ?/max, or equivalently at large 
Xl. We should still have 



drihh 



dr] 



»?=J/max+r)' 



d'pTf{{pT/mN)e'^',PT)^fuW) (3) 



where the universal function fu{v') is independent of ^/s (modulo the separating scale de- 
pendence). So far this is what the experimental data seem to show in both hadron-hadron 
collisions and the heavy-ion collisions. 

Physically, the existence of the limiting fragmentation is a consequence of having a uni- 
versal large x distributions in high energy hadrons combined with the short interaction range 
in the rapidity space Therefore, learning about the limiting fragmentation is equivalent 
to learning about the universal large x distribution. 

In the popular Venugopalan-McLerran model of gluon dynamics, these large x partons 
then act as the color source that generates the small x partons. Therefore establishing the 
validity and also the form of the limiting fragmentation in heavy ion collisions can provide 
an important input for the bulk dynamics of the soft degrees of freedom. 

As far as we can determine, the second universal curve in the transition region has never 
been studied before. In the following sections, we will argue that the appearance of the 
universal transition curve may be anticipated. However, further study is needed to uncover 
the true cause for this universality. 

In this context, it is quite interesting that the deuteron-gold (d+Au) result contains the 
same fragmentation and transition region curve as the gold-gold (Au+Au) result. This is 
discussed in more detail in section IllBI 

The rest of this paper is organized as follows. In section |Hj we analyze available RHIC 
data. A simple parametrization of d'^n/dTf is presented and its consequences explicitly 



calculated. The results from several theoretical models including HIJING|20|, UrQMD 
I22II and a saturation model 
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23L |2J] are compared against the universal curves. Using the 
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FIG. 1: dPn/dri'"^ numerically calculated from the PHOBOS most central collision data[l|. The 
triangles are for ^/s = 200 GeV, the squares are for ^/s = 130 GeV and the circles are for ^/s = 
19.6 GeV. Also shown are two choices of limiting fragmentation functions as explained in the text. 
The two arrows mark the starting point of the plateau {r]p — ymax) and the starting point of the 
fragmentation region (r// — i/max) for y/s = 200 GeV curve. 

two parametrizations of dn/drj from previous sections, we make a prediction for LHC in 
section IIIII Discussions and Conclusions are given in section IIVI Appendix A contains 
details of a calculation not shown in the main text. In Appendix B, we discuss the validity 
(or the lack of) the Wood-Saxon form of dn/drj sometimes used to describe the data. 

II. EXPERIMENTAL LIMITING FRAGMENTATION AND TRANSITION 
CURVES 

A. Analysis of RHIC Au+Au 

If the universal behavior indeed extends up to 90% of the plateau height P], the height 
as well as the total multiplicity would be largely determined by the limiting function fuiv') 
where r]' = r] — ?/max- In reality, the fragmentation region extend up to 50 % of the plateau 
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FIG. 2: (Pn/dr]"'^ numerically calculated from the PHOBOS most central data 'ij. Here rj" = rj — rjp 
where r]p is the location of the hump in dn/dr]. The solid line is guW) = — 0.65r/". 

height at RHIC energy. This fact is hard to see when comparing dn/dri's but becomes 
apparent when comparing d^n/drf^s. In Fig.Q we plot d^n/drf for the most central collisions 
for = 200 GeV, a/s = 130 GeV and = 19.6 GeV numerically calculated from the 
PHOBOS data^. One can see that there are three distinct regions (we will ignore the hump). 
The limiting fragmentation region lies to the right of the minimum of d'^n/drf^ [rj > rjf) in 
which all data points merge together. To its left comes the transition region between the 
fragmentation and the plateau {rjp < rj < r]f). The zero of d^n/dyf is where the plateau 
begins [rj = rjp). This is also the location of the hump in dn/drj. 

It is clear from this figure that the true limiting fragmentation region starts from about 
half way between the plateau and Umax- The area of the triangular shape is the height of the 
plateau. Therefore at these energies the limiting fragmentation region extends up to about 
50 % of the maximum height. Apparent matching of data points below rjf seen in dn/dr]' is 

^ It is not possible to estimate experimental error bars for the slope without knowing the correlation between 
the errors. It is likely that the errors in the neighboring bins are highly correlated. In this paper, we 
assume that this is the case. 
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FIG. 3: Snjdrl'^ numerically calculated from the PHOBOS most central data and semi-central 
data J,]. Here r/" = tj — rjp where rjp is the location of the hump in dn/drj. The solid line is 
gc/(7?") = -0.65r?". 

due to the slow change in the slope but it is not a true universal behavior. 

What is even more interesting is that the transition region also exhibits a universal 
behavior. This is easily seen if one matches the zeros of (Pn/d7f curves (locations of the 
hump in dn/drj) as shown in Fig. |21 One can see that all data points again merge together. 
We will denote this 'universal curve' as 



9u{ti") = 



d^n 



dTf 



{Vp<V< Vf) (4) 

V=Vp+V" 



In Fig. El we also show the semi-central data from PHOBOS together with the central 
collision data. The quality of the data is not as good as the central collision data, but the 
universal behavior is still evident. We do not plot very peripheral data in Fig. |H1 since the 
participant scaling seems not to have been well established for them 0]. Instead, in Fig. U 
we plot the result of pp collisions at various high energies as measured at CERN together 



with a UrQMD calculation from Ref. [25[ and HIJING results. The quality of data for these 
measurements are not as clean as RHIC data from PHOBOS. However, there is a strong 
indication that there is a common transition curve. There is also an indication that the slope 
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□ CERN: ppbar@ 200 GeV 
OCERN: ppbar @ 546 GeV 
ACERN: ppbar @ 900 GeV 
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► HIJING: ppbar @ 200 GeV 
*HIJING: ppbar @ 900 GeV 
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FIG. 4: d?n/dri' numerically calculated from CERN data as compiled by Particle Data Group 2£|. 



Here r/" = rj — ijp where r]p is the location of the hump in dn/dr]. The solid line is guiv") = — 0.40r/". 



in pp (—0.4) is different from the heavy ion result (—0.65). As argued in Section FlI B| this 
is most hkely due to the nuclear modification which is also supported by the HIJING and 
UrQMD results shown in Fig. HJ From now on, we will focus our attention on the central 
heavy ion collisions. 

The shape of dn/dt] is determined by the functional forms of f'uif]') = dfu/drj' and guiv") 
and the condition that these two curves meet at the transition point rj = rjf-. 

guiVf - Vp) = /M-l/max + Vf) (5) 

This is the condition that connects the behavior of the fragmentation region to the plateau 
region. Once the value of rjp is determined by the zero of d'^n/dtf, the pseudo-rapidity 
distribution dn/dr] is fully determined by fu, gu and the condition (0). 

A question then arises: What are the functional forms of the limiting fragmentation curve 
fu and the transition curve gul For the transition curve gu, the current RHIC data shown 
in Fig. 121 suggest that it is a linear function of rj with a y/s independent slope. In this paper, 
we take this to be true and write 

guiv") = -Ki' (6) 
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FIG. 5: The tail part of the PHOBOS data on dn/dr]' for the central 6% of Au+Au collisions at 
^/s = 19.6 GeV. The straight line is {dn/drj') = 0.492 e-^-^^aV. 

where rj" = t] — rjp. The value of K we use is set to 0.65 which is the slope of the straight 
line shown in Fig. (P)). 

As for the limiting fragmentation function, there is little doubt that fu is exponential for 
V > ?/max_as can be clearly seen in Fig. El But what about below ?/max? A current theoretical 
analysis I ISjl relates fu to the gluon distribution function at large x. At moderate Q^, the 
gluon distribution function has the form 

xG'(x,Q2) (7) 

With X = e"^"^™''""'"^ = e~^^'''Q], this means that the limiting fragmentation function 
should behave like 

Mv')-^e-'^'{l-e-'+^'r (8) 

Here —2 in the exponent is due to the mass difference between a proton and a pion. The 
behavior of the expression (jSI) is different from the exponential behavior shown in Fig. in 
the T]' > (rj > ?/max) region. However in the 77' < (r^ < ?/max) region, 

fuiv') ~ e-^"' (9) 
gives a reasonable description with A ~ 0.25jl^. 
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FIG. 6: Pseudo-rapidity distribution for 6% most central Au-Au collisions at ^/s = 200 GeV. Data 
are from Ref. P|. Lower curves are cPn/drf calculated numerically from dn/dr] data. The solid 
lines are our fits. Here we used 9p{x — xq) = 1/(1 + e~^^~^'^^^P) with p = 0.25 and the parameters 
are set to p = 0.95, 1/r = 0.308, q = r /5, 5 = 0.3 and K = 0.65. The value of p is different 
from 1.08 quoted in the text because a finite value of p compensates it a little. Here r]p is not free 
but fixed by the position of the hump. They are at rj' = —3.96, rj' = —3.65 and r]' = —2.6 for 
^/s = 200, 130, 19.6 GeV respectively. 

Combined, the above analysis indicate that the behavior of dn/dr] changes from one 
exponential form to another exponential form when 7] crosses ?/max (or rj' crosses 0) . We may 
represent such behavior with 

r, I dfu P 

^""^"^ ' - df}' ~ eWr + eW? ^^^^ 

In Fig. ^ the solid curve corresponds to this form fitted to rj' > —1 portion of the combined 
data set. The dashed curve corresponds to the extension of the exponential from 77' > 
region. 

If the data points shown in Fig. follow the true universal curve, then we have no choice 
but to conclude that /[/ changes its behavior once rj = z/max is crossed. On the other hand, 
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FIG. 7: dn/dr] and d'^n/drf calculated using Kharzeev and Levin's result . Open symbols are 
PHOBOS results and full symbols are the calculations. 

we may also consider that RHIC energy is not high enough for the true universal curve to 
manifest and what we see in the current data is an accident. One such example is presented 
in Appendix ^ For reasons explained in the Appendix, this accident is unlikely. However, 
only further experiments can give a definite verdict. 
Using Eq. (jH)) for gu^ we write for > 



P 



e,{v' + xf) - K{v' + xp) Op{-Xf - v') OpW + Xp-S) (ii) 



QV'/r _)_ QV'/g 

and 6n(x) is the smeared 6 function with lim 



Opix) = e{x). The 



where r]' = r] - y^ax "x.v. 
minimum of (fn/dr]'^ is located at t]' = —Xf = Vf ~ Umax and the hump of dn/dr]' is located 
at T]' = —Xp = f]p ~ Viaax- The parameters 6 and p control the height and the width of 
the hump. The parameter K = 0.65 is the slope of the transition region. By fitting the 
combined rj' > —1 data, we get p = 1.08, 1/r = 0.308 and 1/q = 1.566 ~ 5/r. We mention 
here that this 1/r value is fairly close to the value of a similar coefficient obtained from 
saturation model studies 
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2J,|27|,|28, 



29[. The shape of dn/dr] obtained from Eq. ()ll|) 



as well as d'^n/drj^ itself is shown in Fig. El In the figure, p = 0.95 is used because having 
finite p changes the slope a little. 
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FIG. 8: d'^n/drf' calculated using Kharzeev and Levin's result|2J| up to the LHC energy. 

The transition between the fragmentation region and the transition region happens at 
7]' = —Xf where the two Hnes meet. This yields the following condition 

K(x.-Xf)- ^.,„/^^,„„ ''pe»'^ (12) 

where the approximation works for large X//'"- This is the condition which relates the 
limiting fragmentation to the plateau and ultimately determines the size of the fragmentation 
region. For large Xf/'^^ the solution is given by 

Xf-Xp-r W{e>'^'''p/Kr) (13) 

where the Lambert function solves w = W{w) exp{W (w)) . With the values of the parameters 
from above, the approximation ()13|1 is good within 1 % for RHIC and LHC but not for SPS. 
For future reference, we note that for large w, W{w) ^Inw — Inlnw. Hence 

X/~lnXp<Xp (14) 

Integrating Eq.flllj) from oo to t]' gives the rapidity distribution dn/drj'. Numerical inte- 
gration yields excellent description of the existing data as shown in Fig. IHl Unfortunately, 
the form of /[/(z/') in Eq.()lU|) does not allow analytic integration in general. However, note 
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FIG. 9: dn/dri' and (fn/dr}'"^ from the central PHOBOS data compared with HIJING vl.37 (filled 
symbols) calculations with dE/dx = — 2GeV/fm js^]- 

that 1/q ^ 5/r. Ifl/g = 5/r, the necessary integration can be carried out in the sharp 
^-function hmit {6p 6). The resulting form is analytic but not very illuminating. Details 
can be found in Appendix 1X1 

Now consider the height of the plateau, {dn/difjQ. Note that the fragmentation region 
(Pn/drf behaves exponentially while the transition region d'^n/dTf behaves linearly in rj. 
Therefore in the large |/max hmit, the contribution from the transition region dominates in 
{dn/dri)o. Physically, this is what one would expect. At high enough energies, the dynamics 
of the central plateau region and the dynamics of the limiting fragmentation region should 
decouple and the height of the central plateau should not depend much on the exact form 
of fu- The height of the plateau in the large j/max hmit is then given by 

This implies 
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FIG. 10: (Pn/dr]''^ numerically calculated from the central PHOBOS data compared with HIJING 
vl.37 calculations. 



since rjp < r]f < ?/max and Xf 
obtained as 



InXp- Integrating once more, the total multiplicity can be 



K 



^total 



maxy 



which implies 



2^1 3, r/ ^ 



(17) 



(18) 



2 

max 



From Eqs. p5j) and (fTTjl we conclude that the central plateau cannot rise faster than y 
or In^ s and the total multiplicity cannot rise faster than y^^x ^- The only possible 

way to get faster dependence is to have a/s dependent K, or faster rising gu (for instance 
an exponential). Judging from Fig. HI this is not likely up to 900 GeV. Also there is an 
additional evidence from the CDF collaboration |^] that up to i/i = 1.8 TeV, the central 
plateau in pp collisions rises only as fast as In^ s. 

In many current models of heavy ion collisions, ni,nta,i grows faster than In^ s. For instance, 
Ref. [2^ has ntotai ~ s^^'^ and the e~^e~ model[l3| has ntotai ~ gcvlns -^i^QYe X and c are 
constants. Parametrization of pp, pA and AA data up to the SPS energy by Gazdzicki 
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FIG. 11: d^n/dr]''^ numerically calculated from the central PHOBOS data compared with UrQMD 
calculations (filled symbols). The solid line is UrQMD result for y/s = 200 GeV without rescatter- 
ings. 

and Hansen gives Net oc s^f^. At present energies, these are indistinguishable from 
polynomials in Ins. However, as will be presented shortly, LHC will be able to tell whether 
the bound (fTHj) indeed holds for high energy heavy ion collisions. 

At this point, we can attempt a partial explanation of the appearance of the universal 
transition curve. Suppose that as the collision energy becomes larger the dynamics of plateau 
region largely decouples from the dynamics of the fragmentation region. This is certainly 
the case for the fu and gu given in this section as indicated by Eq. lfT^ . Eq. lfT^ implies 
that in the large j/max limit, Tjj ^ rjp and 

??/ = Z/max + O(lnymax) • (19) 

This is a consequence of having an exponential fragmentation curve and a polynomial tran- 
sition curve. Note that the functional form of gu and fu enters only through the logarithmic 
correction. 

At RHIC energies, the area under fu and gu looks like an isosceles triangle. This is 
because Inj/max is still not that small compared to ?/max- However since an exponential rises 
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FIG. 12: The derivative {(fN/drf)/{Npi,^t/'i) with respect to ry. Data are taken from the PHOBOS 
website kai. 



fast, the area will look more and more like a right triangle as the energy grows and the area 
will become dominated by the transition part: 




rip 



Vf 



dV9u{v-Vp) 



dv" 9uiv") 



^max 



dv" 9uiv") 



(20) 



/o ^0 

Therefore, to leading order in ?/max! {dn/dri)^ is a function of i/raax ~ ^^{y/s/rriiy) and it is 
independent of the functional form of the fragmentation curve /[/. Denoting the functional 
dependence as {dn/dri)^ = S{ymax), the universality of gu follows if the following relationship 
holds 



guiv") 



dS{i') 
dr]" 



(21) 



Once the dependence of {dn/dr])^ on y/s is given, gu is totally determined and it is indeed 
universal up to logarithmic corrections. The relation ()21|) certainly holds for Eqs.® and 
(Uni) when ?/max > 1- 
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FIG. 13: The same as Fig ll2l but d+Au data is vertically scaled by a factor of 1.3 and shifted to 
the right by 0.4 unit of pseudorapidity (2 experimental bins). 



The hole in this argument is that the relationship pO|) does not automatically imply 
Eq.(PT|). For instance, suppose S{yma.x) = c^Z/max- this case, any 



guiv') 



^ l/max 



in 



//\n-l 



(22) 



with n > 1 satisfies Eq. ipnj) . Unless n = 2, however, guiji") depends on ?/max and hence it 
is not universal. Surprising fact is that the data seems to suggest n is indeed 2 or at least 
very close to it. 

The relationship (j2H) is remarkable. It relates an observable that is a function of collid- 
ing energy to an observable that is a function of the pseudo-rapidity at any fixed energy. 
Unfortunately, energies probed so far are too small for this to manifest. As seen in Figures 
121111 the transition region is not truly dominant yet. However, we should be able to test this 
relationship at LHC 

It is also instructive to compare some theory curves with RHIC data as shown in Figs.O- 
im As shown in FigslT] the saturation model by Kharzeev and Levin 0] gives a good 
description of the plateau and the transition region at RHIC energy although the fragmen- 
tation region is badly off. However, since the model is based on small x picture, it is not 
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FIG. 14: The same as FigUHlbut for dN/drj. 

supposed to be valid in the fragmentation region. From the expression of dn/dy given in 
Ref . [2^ , it is clear that the transition curve obtained by Kharzeev and Levin is exponential 
and this form of dn/dy does not satisfy the relation ()2ip. The reason Fig. |H1 shows approxi- 
mate universal behavior up to ^/s = 1000 GeV is A is small. At LHC energy, the violation 
of the universality is clearly seen for this model. 

HIJING results js^js^ with shadowing and a parton energy loss of dE/dx = —2 GeV/fm 
and an energy dependent scale parameter po as considered in [3^ are shown in FigEl It 
is quite evident in the d'^n/drf plot that the fragmentation region dominates in HIJING. 
Again to test the transition curve universality, one must go beyond the RHIC energy. FiglTUl 
shows d'^n/dTf up to the LHC energy. From the figure, it is quite clear that HIJING 
does not contain a universal transition curve. Furthermore, at higher energies, the central 
region develops a bump instead of a plateau. This feature is due to the abundance of the 
minijets. As can be seen in Fig. El by enhancing the parameter pq ( equivalent ly, reducing 
the number of minijets) HIJING becomes closer to the other models. But the transition 
region universality is clearly not a feature in the HIJING model. 

On the other hand, it is quite striking that the default UrQMD results get both dn/drj and 
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FIG. 15: The same as Fig ll2l but d+Au data is vertically scaled by a factor of 1.9 and shifted to 
the right by 0.2 unit of pseudorapidity (1 experimental bin). 

(Pn/drf right at RHIC energies. It is also significant that without rescatterings, UrQMD 
does not describe the data well. Why UrQMD results in universal transition curve is not 
yet clear. 

B. Analysis of RHIC d+Au 

Recently the PHOBOS collaboration published the result of measuring the pseudora- 
pidity distribution of produced particles in the deuteron-gold (d+Au) coUisions at y/s = 
200 Gevjs^]. At a first glance, it would seem that there is no common feature at all between 
the d+Au dn/di] and Au+Au dn/dt], especially if one just looks at the participant-scaled 
results. However, when dealing with very asymmetric systems such as d+Au, one must be 
careful about the scaling behavior. As can be easily shown in a simple wounded nucleon 
model, the scaling of produced particles in the heavy ion side and the d side should be 
different. The number of wounded nucleons in the heavy ion side depends on the linear 
size of the heavy ion whereas the d side always have 1 or 2 wounded nucleons. Hence, the 
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FIG. 16: The same as Fig ll5l but for dN/dr]. d+Au result is shifted vertically down by 3.6. 

multiplicity in the heavy ion side should have an additional scale factor ~ A^^^ compared 
to the d side. 

To see whether there is a common feature between Au+Au and d+Au or not, again it is 
much better to look at the derivative {(f N / dr]"^) / {Np^rt / 2) as shown in in FigHSl Judging 
from this figure, it is clear that there is a common feature. To bring it out more clearly, we 
vertically scale the d+Au result by a factor of 1.3 and shift it horizontally by 0.4 unit of 
rapidity (or 2 experimental bins). This results in Figs. El and El which leaves no room for 
doubt that the shape of dN/dr] for r] > 1.5 is common to both Au+Au and d+Au results. 
It is also interesting to see that different scaling (additional factor of 1.5 compared to FiglT^ 
and the rapidity shift of 0.2 (or 1 experimental bin) instead of 0.4) brings the Au side of 
the spectrum together as shown in FiglT^ and FiglTHl Again, there is no room for doubt 
that there is a common curve. This implies that beside a constant component, the shape of 
dn/dr] for both Au+Au and d+Au is simply related by scaling even in the Au side. 
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FIG. 17: Predictions for central collisions at LHC. Also shown are parametrized RHIC 200 result. 
We set y/s = 5.5 TeV and Pb+Pb for LHC. 

III. PREDICTION FOR LHC 

Given the two forms of parametrization considered in the previous sections, it is possible 
to extrapolate and predict what should happen at LHC where ^/s = 5500 GeV and z/max = 
8.68. To do so, we need to parametrize the functional form of rjp or rjf. The value of rjp in the 
model from section 111 Al is not a free parameter. The position of the hump clearly visible in 
dPn/drf as a zero is the value of r/p. From the PHOBOS data, one gets, ymax^Vp = 3.96, 3.65 
for ^/s = 200, 130 GeV. For i/i = 19.6 GeV case, the data has ?/max — Vp = 2.75. However, 
with our model, ymax — Vp = 2.6 describes the data better. 

By fitting the above values of rjp with ?/max — Vp = ^U'^ + p\ny and ?/max — Vp = "^U^ + C*, 
we get^ 

^ f 0.60 + 0.73^^91 (Model I) 

^"'^ 1 0.33 In + 0.96 y^I^ (Model H) ^ ' 



Since we have only three data points, one cannot fit the full ^max ~ Vp = + /31ny + C. 
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2/max Vp 


{dn/dT])o 


'^total 


Model I 


5.8 


6.9 


87 


Model II 


5.6 


6.5 


83 


K & L 




10.7 


110 


HIJING w/ po = 3.5 GeV/c 




21.4 


160 


HIJING w/ po = 7.0 GeV/c 




11.6 


100 



TABLE I: Predictions for LHC central collisions. We set ^ = 5.5 TeV and Pb+Pb for LHC. 



These two parametrizations do not differ much up to ?/max = 10. The LHC predictions from 
these two parametrizations are given in Table U. The shape of dn/dt] is given in the FiglTTI 
together with the Kharzeev- Levin prediction and the HIJING predictions with two different 
minimum minijet energies. The results obtained for pq = 3.5 GeV suggested in reference 
js^ . is clearly very different from other models. Increasing the mini-jet scale parameter to 
a higher value po = 7.0 GeV brings it to a better agreement with other models. Only data 
from LHC will allow us to draw a definite conclusion and to choose the right value for this 
parameter p^. 

One striking feature is that in both the Kharzeev- Levin model and the HIJING model 
with pq = 3.5 GeV, the central plateau disappear. This is mainly due to the fact that 
these models do not contain rescatterings of the secondaries and hence cannot not undergo 
a Bjorken-like expansion. 



IV. DISCUSSIONS AND CONCLUSIONS 



In this paper, we showed that there exists a second universal behavior in the rapidity 
distribution of produced particles. The data we have analyzed clearly indicate that d^n/drf^s 
in the transition region taken at different energies follow a common curve. This is not easy 
to see when comparing dn/drj^s but clearly seen when comparing d^n/diff^s. 

We emphasize here that any model that purports to describe the rapidity distribution 
in the whole rapidity space must be able to reproduce not only the limiting fragmentation 
curve, but also the universal transition curve. 

The existence of the two universal curves implies that the shape and the size of the 
rapidity distribution itself is mostly determined by (i) the limiting fragmentation curve fjj, 

^ To predict the height more accurately, we need to remember that Ea. (|15|l neglects some part of the tail 
contribution. 
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(ii) the universal transition curve gu and (iii) the starting point of the plateau region rjp. 
Non-trivial physics resides in the y/s dependence of r]p or equivalently the size the central 
plateau. 

The physics behind the limiting fragmentation curve is well known to be the Feynman- 
Yang scaling which states that at high energy, large xl behavior of an inclusive cross- 
section is independent of ^/s. In this paper we argued that the physics behind the universal 
transition curve is in fact the decoupling between the dynamics of the plateau and the 
fragmentation region. We note that it is also intriguing that perhaps a connection to the 
gluon parton distribution can be made. In saturation models, dn/drj is related to the gluon 



distribution function 



18l I23I |24| . In this case, the universal transition curve puts a severe 



restriction on the behavior of the gluon distribution at moderate x. 

In this work, we found that the universal transition curve is linear in t]" = rj — rjp based 
on -^i = 20 — 200 GeV RHIC data and UA5 pp data. A consequence of having a linear 
transition curve is that the plateau height {dn/dr])o cannot grow faster than In^ s and the 
total charged multiplicity cannot grow faster than In^ s. This polynomial behavior in Ins is 
maintained if the transition curve is polynomial in 77". A power law growth {dn/dri)o ~ 
or an e~^e~ type exponential growth is possible only if the transition curve is exponential. 
The available data does not show such exponential behavior. It does not, of course, rule out 
a change in the behavior at higher energies. As shown in the last section, these possibilities 
can be clearly distinguished at LHC. 

What we have found in this study also impacts hydrodynamic studies. As can be seen 
in Eq.(j2SI) the size of the plateau does not grow fast. Moreover as Umax 00, rjp/y^i^^ — > 0. 
Therefore, the relative region of validity for 2-D hydrodynamic calculation shrinks as the 
energy goes up and the need for 3-D hydrodynamic calculation becomes greater. Further- 
more, the existence of the universal transition curve will tightly constrain the longitudinal 
evolution of the hydrodynamic system. 

We have also analyzed the deuteron-gold result from RHIC and found that there is a single 
common curve that determines the shape of dN/ dt] for both d-l-Au and Au-|-Au cases. A few 
conclusions can be drawn from our analysis. First of all, the different scaling factors for the 
deuteron side and the gold side indicate that the scaling of d-l-Au system is more complex 
than a simple participant scaling. This implies that using a simple participant scaling can 
potentially mislead the comparision between the d-l-Au result and Au-|-Au result. This is 
especially significant for the Au side where there appears to be a constant component on 
top of the two universal curves discussed in this paper. Again, these facts are much more 
transparent if one compares d'^N/dTf. 

Second, the existence of a function common to both Au-|-Au and d-l-Au indicates that the 
dynamics of the transition region and the fragmentation region in the Au-|-Au case cannot 
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depend much on the final state interactions. It can only depend on initial state parton- 
nucleus dynamics. Especially, whether or not a hot and dense system is formed in Au+Au 
collisions does not infiuence the shape of dN / dr] beyond the plateau region. 

The exact physical meaning of the rapidity shifts and the constant component in the 
d+Au data are under investigation. 
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APPENDIX A: INTEGRAL OVER dfu/drj' 

To calculate dn/dr], we need 

£ T^ = -P £ eW. + e5.'A (Al) 
Changing variable to y = e^" 1'^ yields 
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These integrals can be found in integral tables, for instance Ref. 42|. We get 

// dfu J -n'/r , P 



dr]" 



dr]" 



pr 



e-'^'l' + ^ (tan-i (l + v^e"'/") - tan" 
_ilo / l + v^e"'''^'' + e'''''^'^ \ 1 



1 - v^e"'/' 



TT 



(A3) 



APPENDIX B: UNIVERSAL FRAGMENTATION CONDITION FOR WOOD- 
SAXON 

A popular choice of parametrization for dn/drj is the Wood-Saxon function. Many models 



for hadron-hadron collisions developed in the 70 's [36|, |37|, |38|, [s^ also had this type of 
dn/drj. As will be shortly shown, this Wood-Saxon form is not compatible with the transition 
region universality and hence it is unlikely that this is the right form of dn/drj. Nevertheless 
we feel that it is worth considering the Wood-Saxon form here because it gives an example 
of slowly changing (as opposed to universal) limiting fragmentation curve. 

A reasonable description of the current data can be provided by the following combination 
of the Wood-Saxon (Fermi-Dirac) functions and a hyperbolic cosine 41 1 



dn 
drj 



g cosh(?7/C) 



(Bl) 



where rjj, ( and cr are functions of ?/max- Here the parameter rjj roughly corresponds to 
where the fragmentation region begins. The hyperbolic cosine is there to provide the dip in 
the middle. Since the dip is usually shallow, C ^ cr- 

Since we are not so much interested in the dip, we consider a simplified form"^ 



dn 
dr] 
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for ?7 > 



1 + fAv-Vf)/^ 

Universal fragmentation behavior demands that for high enough energy 



(B2) 



dn 
drj 



fuW) 



(B3) 



'7=ymax+'?' 



where fu{v') is indepdent of ?/max- i]f /a ^ 1, this just implies that so 
that it compensates the large exponential in the denominator. Near rj = ?/max this yields 
f^{j^') = Koe-'''/". However at SPS and RHIC energy, y^a.x is only about 3 to 5 and it can 
be easily shown that simply having g = /tge^^™'"''"''^'^^'^ with a constant a does not result in 



Entirely analogous anaylsis can be also performed using Eq.(jBTj). 
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the universal limiting fragmentation. Instead, we must regard all parameters appearing in 
Eq. ()B2|) as functions of ?/max and look for a relationship among them by requiring 

C'Z/max dr]' 

near r] = ?/max (or 7]' = 7] - ?/max = 0). 

The solution of Eq. ()B4|) is obtained as follows. Set a = 1/t and rj = y^ax + v' 
Xf = Z/max -Vf and write 

rdn_ g 

^ - dr]' 1 + e('''+^/)- ^ ^ 

Taking the derivative with respect to y = Umax yields 

a/ ^ {dg/dy){l + e(^'+>^/)n - ^ [v'idr/dy) + ^(x/t)/^!/] e(^'+>^/)- 

(1 + e(v'+xf)ry ^^'^^ 

If this is to be independent of r]' for small rj', we must have 

(dg/dy) [1 + e-^f^{l - r^'r)] ^ g [v'idr/dy) + d{xfr)/dy] (B7) 
which yields the following two conditions: 

^dq dr ,^ , 

ay ay 



and 



Assuming monotonic functions, we can rewrite them as 

dg dr 



h + e--r] = g^i^ (B9) 



e^f^ (BIO) 
9 T 



and 

c^fl' _ d{xfr) 



g 1 + e"^/" 
Solving the second equation first gives 



(Bll) 



g = Ko{l + e^?'— -''/)/'^) (B12) 
Combine the two equations to get 

(1 + e-^/^) T ^ ' 
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Let L = e and use dL = —dixfT) e ^^'^ to get 

dL dr 
~ l+L ~ ~~ 

Solving this equation, we finally get 



(B14) 



r = C"(l + e-'^^") (B15) 

which can be rearranged to yield 

a = ^ ^ (B16) 

or 

l = l + ^W ((x//^o)e-^^/'^«) (B17) 

where we used X/ = 1/max ~ ''7/ and W{vS) is the Lambert function that solves w = xe^ 
for X. Hence given a value of r]f, the width and the height of the Wood-Saxon function is 
completely determined. 

For large enough (y^ax - Vf)/^^, 



1 1 

_ ~ —(1 + e~*^^™^''~''-f)/°'° 
a (To 

ignoring terms of 0(e~^*^^™'"''~^^-'/°""). At an asymptotically high energy 



) (B18) 

a (To 



lim a = (To (B19) 

Hence the limiting curve is given by 

fu{v') = ^oe-'^'/'"' (B20) 

Fitting the r] > ?/max portion of = 19.6 GeV data yields 1/ao = 1.253 and kq = 0.492 
as shown in Fig. El The height of the plateau and the total multiplicity can be now easily 
obtained from Eq. ()B2|) 

/dn\ ^ /to(l + efa— ''z)/-) 

\dri)^ 1 + e-''//'^ ^ ' 

and 

ntot.i = 2r^f (^^^ +Oie-^f/n (B22) 

The resulting dn/dr] and d'^n/drf^ are shown in Fig. IT^ together with PHOBOS data. 
Note that although fragmentation region universality is reasonably well described by the 
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FIG. 18: {Xf,f^} = {1.32,0.694} for 200 GeV, {1.20,0.681} for 130 GeV and {0.77,0.619} for 
19.6 GeV. The limiting curve parameters are kq = 0.492 and l/ao = 1.253. 

Wood-Saxon functions, the transition region universality is not. Again, we emphasize that 
it is the slope {(Pn/drf) that gives clearer criterion for the goodness of the description. 

Unlike the previous case, the Wood-Saxon case has no separate rip. This is because 
both the fragmentation and the transition behaves like an exponential ~ e'~\'^~'^f\/'^° near 
r]f. Therefore the transition between the plateau and the limiting behavior happens within 
about 3a"o around 77/. This fact also indicates that dynamics of the plateau and the dynamics 
of the fragmenation region does not decouple even at an asymptotically high ^/s. 

This non-decoupling also allows us to put a severe condition on the transverse energy. 
To calculate the energy content of the plateau, we need to carry out an integral over the 
product of the Wood-Saxon and a hyperbolic cosine. This can be done, but the resulting 
form is not particularly illuminating. However, within the plateau we can approximate 

( dn\ 

^plateau ~ 2 j {mr)^,, sinh(r7/) (B23) 



Energy conservation demands that 



(^^) < e-(V^-i)(2/ma.-^/) (B24) 
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ymax - Vf 



{dn/dr])o 



"-total 



WS I 



2.3 



10.6 



130 



WS II 



2.1 



8.5 



110 



TABLE II: Predictions for LHC central collisions using the Wood-Saxon form. We set y/s = 5.5 TeV 
and Pb+Pb. 

where we used Eg. ()B2Hl . Since a < 1, this indicates that the average transverse energy in 
the plateau region must be a decreasing function of Umax if Vmax — Vf is an increasing function 
of ?/max- This is an absurd resuh. One would expect that as becomes larger, (^t)pi 
would also become larger or at least reach a limiting value, but not decrease. This, in our 
opinion, invalidates the Wood-Saxon description of dn/drj. 

Nevertheless, it is instructive to also have the extrapolated Wood-Saxon result to LHC. 
For the Wood-Saxon form, we find that Xf = 1-32, 1.20,0.77 for ^ = 200 GeV, 130 GeV, 
20 GeV, respectively. These yield 



The results for LHC are tabulated in Table ITTl The values of {dn/difiQ and ntotai are compa- 
rable to the saturation model (K & L) values in Table lU Also one can easily see in Fig. 
that the limiting fragmentation curve followed by the Wood-Saxon functions is not the same 
as the one followed by the interpolating-exponential ones. 




(WS I) 
(WS II) 



(B25) 
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FIG. 19: Predictions for central collisions at LHC. Also shown are parametrized RHIC 200 result. 
We set ^/s = 5.5 TeV and Pb+Pb for LHC. 
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